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A phase transition from paramagnetism to ferromagnetism in neutron star interior is explored. 
Since there is 3 P2 neutron superfluid in neutron star interior, it can be treated as a system of mag- 
netic dipoles. Under the presence of background magnetic field, the magnetic dipoles tend to 
align in the same direction. Below a critical temperature, there is a phase transition from para- 
magnetism to ferromagnetism. And this gives a convenient explanation of the strong magnetic 
field of magnetars. In our point of view, there is an upper limit for the magnetic field strength of 
magnetars. The maximum field strength of magnetars is about (3.0 — 4.0) x 10 15 G. This can be 
tested directly by further investigations. 

Magnetars are instable due to the ultra high Fermi energy of electrons. The Landau column 
becomes a very long cylinder along the magnetic field, but it is very narrow and the Fermi energy 
of electron gas is given as E F (e) &A0{B/B cr ) l l A when B>B„. E F (e) » 90MeV Whenfi ~ 10 15 
G. Hence, the electron capture process e~ + p — > n + V e will be happen rapidly. Thus the 3 P2 
Cooper pairs will be destroyed quickly by the outgoing neutrons with high energy. It will cause 
the isotropic superfluid disappear and then the magnetic field induced by the 3 P2 Cooper pairs 
will be also disappear. These energy will immediately be transmitted into thermal energy and 
then transformed into the radiation energy with X-ray - soft y-ray. We may get a conclusion that 
the activity of magnetars originates from instability caused by the high Fermi energy of electrons 
in extra strong magnetic field. 
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1. Introduction 

A puzzle is what is the origin of strong magnetic field of magnetars and it is still open question, 
although it has been investigated by many authors. 

In this paper we propose a new idea for the origin of the magnetars. The strong magnetic fields 
of the magnetars may originate from the induced magnetic moment of the 3 P2 neutron Cooper pairs 
in the anisotropic neutron superfluid. 



2. Induced paramagnetic moment of the 3 P2 neutron superfluid and the Upper limit 
of the magnetic field of magnetars 

A magnetic moment of the 3 P2 neutron Cooper pair is twice that of the abnormal magnetic mo- 
ment of a neutron,2/x„ in magnitude, and its projection on the external magnetic field (z-direction) 
is a x (2jj. n ) , a z = 1,0,-1, where \i n is the absolute value of the magnetic moment of a neutron, 
\i n = -0.966 x 10~ 23 erg/G. 

A magnetic dipole tends to align in the direction of the external magnetic field. The 3 P2 neutron 
Cooper pair has energy Oz x pL n B in the applied magnetic field due to the abnormal magnetic 
moment of the neutrons. B is the total magnetic field including the background one. 

The difference of the number density of 3 P2 neutron Cooper pairs with paramagnetic and 
diamagnetic moment is 

An T =n_ l -n +i =n n ( 3 P 2 )f(^) (2.1) 

= 2sinh(2x) 
J v ' l+2cosh(2x) 

The Brillouin function, f{pL n B/kT), is introduced to take into account the effect of thermal motion 
We note that / is an increasing function, in particular,/(x) rj 4x/3, for x -C 1 and f(x) — > 1, when 
i>1. f{pL n B/kT) increase with decreasing temperature. And this is the mathematical formula for 
the B-phase of the 3 P2 superfluid. 

A relevant question is how many neutrons have been combined into the 3 P2 Cooper pairs? 
The total number of neutrons is given by N = 2V Jq f = ^k\ (Here k is a wave vector). 
The neutrons combined into the 3 P2 Cooper pairs are just in a thin layer in the Fermi surface with 
thickness k&^, = \/2m„A, 1 ( 3 / 5 2) (We would like to emphasize that the energy gap, A, is the 
binding energy of the Cooper pair rather than a variation of the Fermi energy due to the variation 
of particle number density). Then we have SN « ^?^fj(^A << kp). Thus, the fraction of the 
neutrons that combined into the 3 P2 Cooper pairs is 



SN k A 

The fraction of the neutrons that combined into the 3 P2 Cooper pairs is (p = Tik) 

_ 47r^[2m n A( 3 P 2 )] 1/2 _ or A( 3 P 2 ) ll /2 n 
q ~ (4n/3)pl ~" [ E F 1 ■ {ZA) 

The Fermi energy of the neutron system is Ep « 60 {-^—) 2 ^ Me v. The energy gap of the 
anisotropic neutron superfluid is A( 3 P 2 ) ~ 0.05 MeV |0|, q ~ 8.7%. Thus, the total number of the 
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3 P2 Cooper pairs is n n (^Px) ~ qNAtn( 3 P 2 )/2 . Therefore, the total difference of the 3 P2 neutron 
Cooper pair number with paramagnetic and diamagnetic moment is 

AA^ = «„( 3 P 2 )/(^) = l -N A mCP 2 )qf(^) (2.5) 

The total induced magnetic moment, of the anisotropic neutron superfluid is 

(3p 2) = 2pLn xANt = i^N A m( 3 P 2 )qf(^B/kT) (2.6) 

Where m( 3 P2) is the mass of the anisotropic neutron superfluid in the neutron star, N A is the Avo- 
gadro constant. The induced magnetic moment is just the fully magnetized quantity pL n NAm( 3 P 2 ), 
with two modification factors. The factor q takes into account the Fermi surface effect. While 
f(jX n B/kT) is the thermal factor taking into consideration the finite temperature effect. For a dipo- 
lar magnetic field [Q] |jUjvs| = ByR 3 NS /2. Here B p is the polar magnetic field strength and R^s is the 
radius of the neutron star. The induced magnetic field is then 



B (in) = Zr_^l = "^^^' qf^B/kT) (2.7) 
%S K NS 



2 ^ in \ 3 P 2 ) = I^NaC Pi) 

Rle Rl 



bW( 3 P 2 ) = 2 ^ N f Pl) q G « 2.02 x 10 14 T] G (2.8) 

R NS 



m?Pi) 3 , A„( 3 P 2 
0AM sm NS > 6l 0.05MeV 



Here tj is the dimensionless factor describes both the macroscopic and microscopic properties 
of neutron stars. 

We note that the induced magnetic field for the anisotropic neutron superfluid increases with 
decreasing temperature as the Brillouin function f(fj. n B/kT) which tends to 1 when the temperature 
decreases low enough. Actually, f{pL n B/kT) ~ 1 as long as \i n BjkT ^> 1. For example, this is true 
when T < 10 7 K if B = 10 15 G. 

There is an upper limit for the induced magnetic field of the 3 P 2 superfluid according to 



eq.([2T8|). It corresponds to the maximum value unity of the temperature factor f(fX n B/kT). This 
upper limit can be realized when all the magnetic moments of the 3 P 2 neutron Cooper pairs are 
arranged with the paramagnetic direction as the temperature become low enough. 

The maximum magnetic field for magnetars depends on the total mass of the anisotropic neu- 
tron superfluid of the neutron star. The upper limit of the mass for the neutron stars is more than 
2M sun [Q] . It is therefore possible that the mass of the anisotropic neutron superfluid of the heav- 
iest neutron star may be about (1 — l.5)M sun . Hence, the maximal magnetic field for the heaviest 
magnetar may be estimated to be (3.0 — 4.0) x 10 15 G. This can be tested directly by magnetar 
observations. 

3. From paramagnetism to ferromagnetism 

We have \i n B < kT when B < 10 13 G and T > 3 x 10 6 K, therefore we could make Taylor 
series 

3 R^ s H kT T 7 ' 
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B = B (in ^+B {0) (3.2) 

Where B^ ' is the applied magnetic field which includes both the initial magnetic field of the col- 
lapsed supernova core and the induced magnetic field produced by the Pauli paramagnetization of 
the highly degenerate relativistic electron gas in the neutron star interior m\. 



We may solve the induced magnetic field by combining eq.(3.1) and eq.(3.2): 



B (m) = rJz_ _ n-iflto) ( 3 3) 

The induced magnetic field, B^ , of the anisotropic neutron superfluid is much weaker than the 
applied magnetic field , B^°\ when the temperature is very high, Tg 3> 1.9T7 . However, the induced 
magnetic field of the anisotropic neutron superfluid would exceed the applied magnetic field when 
the temperature decreases down to 1 < Tj/ (1 .9r) ) < 2. Therefore B^ in ' should be calculated by both 
eq.(2.£) and eq.(3.2) when Tj ~ 1.9x]. This belongs to the domain of ferromagnetism (see [Q][{7|]). 



In view of the above discussions we may get a very important conclusion: the strong magnetic 
field of magnetars may originate from the induced magnetic field by the ferromagnetic moments 
of the 3 P2 Cooper pairs of the anisotropic neutron superfluid at a moderate low temperature about 
10 7 K. 



4. Activity of Magnetars 

The core temperature of the magnetars is about 10 7 K in our model. While observations show 
that some SGR's and AXP's have high thermal-type-spectrum X-ray flux, being among the hottest 
neutron stars. We discuss this question and give reasonable and consistent explanation in this 
paragraph. 

Energy of an electron under the super strong magnetic field is quantized. It is 

E 2 (p z ,B,n,o) = m 2 e c 4 + p 2 z c 2 + (2n + 1 + o)2m e c 2 ]X e B (4.1) 

(The magnetic moment of an electron is 0.927 10~ 20 erg /G) The overwhelming majority of elec- 
trons congregate in the lowest levels n = or n = 1... when B S> B cr . The Landau column is a very 
long cylinder along the magnetic field and it is very narrow. 

We may find the Fermi energy of the electrons by solving the following complicated integral 
equation 



Ntotai = N A Y eP = -3- / d Pz £ / 8{{p x /m e c) - [(2n + 1 + a)b\ l/2 ) P± d P± (4.2) 

n max (p z ,b, a) = lnt{U{-^) 2 - 1 - (-^) 2 ] - 1 - a} (4.3) 
2b m e c z m e c 

b = B/B cr (4.4) 

Ntotai is the total occupied number of electrons in a unit volume and Int{ ... } is the integer function. 
After some calculation, we may get the approximation expression 

E F (e) «40(B/B cr ) 1/4 when (4.5) 
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We have E F (e) « 90MeV, when B ~ 10 15 G. 

Magnetars are instable due to the ultra high Fermi energy of electrons. Hence, the electron 
capture process e~ + p — > n + v e will be happen rapidly as long as the Fermi energy of electrons is 
much greater than 60 Mev which is the Fermi energy of neutrons. Energy of the resulting neutrons 
will be rather high and they will react with the neutrons in the P2 Cooper pairs and thus the P2 
Cooper pairs will be destroyed quickly by the process n+(nf,n|)— >n + n + n. It will cause the 
isotropic superfluid disappear and then the magnetic field induced by the 3 P2 Cooper pairs will be 
also disappear. 

The remaining energy per outgoing neutron after the process above will be approximately 

e « \[E F {e) + E F (p)-E F (n) - (m n - m p - M e )c 2 ] (4.6) 

These energy will immediately be transmitted into thermal energy and then transformed into the 
radiation energy with X-ray - soft y-ray. 

The total released energy may be approximately estimated as 

E^ « £ x n( 3 P 2 ) =3ex qN A m{ 3 P 2 ) » 2.0 x 10 51 [(— ) 1/4 - 1 .5] "f ^ ergs (4.7) 

B cr 0.lM mn 

The x-ray Luminosity of AXPs is L x « 10 34 — 10 36 erg /sec. It will be enough to maintain the 
luminosity of AXPs over 10 8 yr. 

We may get a conclusion that the activity of magnetars originates from instability caused 
by the high Fermi energy of electrons in extra strong magnetic field. However, we have to 
calculate rate of the electron capture process for given magnetic field to see if our idea is reasonable 
for the x-ray luminosity of AXPs and their surface temperature. This is our next work. 
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